Different faces of chaos in FRW models with scalar fields — geometrical point of view 
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FRW cosmologies with conformally coupled scalar fields are investigated in a geometrical way 
by the means of geodesies of the Jacobi metric. In this model of dynamics, trajectories in the 
configuration space are represented by geodesies. Because of the singular nature of the Jacobi 
metric on the boundary set dD of the domain of admissible motion, the geodesies change the cone 
sectors several times (or an infinite number of times) in the neighborhood of the singular set dD. 
We show that this singular set contains interesting information about the dynamical complexity 
\f} • of the model. Firstly, this set can be used as a Poincare surface for construction of Poincare sections, 

and the trajectories then have the recurrence property. We also investigate the distribution of the 
intersection points. Secondly, the full classification of periodic orbits in the configuration space 
£Nj ■ is performed and existence of UPO is demonstrated. Our general conclusion is that, although 

the presented model leads to several complications, like divergence of curvature invariants as a 
measure of sensitive dependence on initial conditions, some global results can be obtained and some 
additional physical insight is gained from using the conformal Jacobi metric. We also study the 
complex behavior of trajectories in terms of symbolic dynamics. 
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I. INTRODUCTION 



In the present paper we work in a certain geometrical model of dynamics which can be constructed due to some 



variational principles of mechanics that allow one to construct an isomorphism between the dynamics and its model. 

Recently a geometric description of chaos in Hamiltonian systems of cosmological origin has been formulated using 
the tools of Riemannian (or pseudo-Riemannian) geometry P, Q • We concentrate on the approach to dynamics of 
^ the systems with a natural Lagrangian function C = \g a fi(l a (f — V(q). On the level of constant energy E this system 
I \ can be reduced to the geodesies flow on a pseudo-Riemannian manifold with the Jacobi metric g a p = 2{E — V)g a p 
[HBHHlii The conformal metric becomes degenerate for certain values of energy E, at some points of configuration 
space {E — V} for classical systems as well as for general relativistic ones. As a consequence one has the complication 
, of matching together geodesies well defined on open sets across a singular surface and the divergence of curvature 
invariants characterizing the property of sensitive dependence on initial conditions. One can avoid this crucial problem 
addressed in the context of Mixmaster models 0, 13 by formulating corresponding dynamical system as a system in 
' the Finsler [j| or Eisenhart metric [Tof . 

The fact that Jacobi geometry is singular suggests that this model of dynamics [34| is one of the worst possible 
choices when it comes to achieving characterization of property of sensitive dependence on initial conditions it terms 
of curvature invariants or when it comes to achieving global results. It will be demonstrated that it is not true for 
the example of the behavior of geodesies in the Jacobi metric for FRW cosmological models with conformally coupled 
scalar field. Of course, because of the singular nature of the Jacobi geometry, the geodesies change from time-like to 
space-like (sometimes several or an infinite number of times) during their evolution, but in order to "piece together 
the results" , one is forced to try to match geodesies (which are representing periodic orbits of their original dynamics) 
each time the solution passes a singular point. This leads to the full classification of periodic orbits in terms of "just 
piecing together" segments of geodesies. It is an example of the fact that some global theorem about geodesies can 
be achieved. Using this classification numerous methods based on symbolic dynamics and detection of the existence 
of unstable periodic orbits are proposed toward an invariant characterization of dynamical complexity in cosmology 



The organization of the paper is the following. In section [H] the cosmological FRW model with conformally coupled 



* Electronic address: 
t Electronic address: 



hrycyna@byk.oa.uj.edu.pl 
uoszydlo@cyf-kr . cdu. pi 



2 



scalar field is presented as a Hamiltonian dynamical system. This system belongs to larger class of simple indefinite 
dynamical systems for which kinetic energy form is quadric in the momenta and has Lorentzian signature (is indefinite) . 
In section ITTT1 we consider Poincare sections of the trajectories, with the surface formed by the boundary set of the 
domain classically admissible for motion, as a indicator of complex behavior. Note that the trajectories have property 
of recurrence due to multiple intersections of singular set. In this section the full classification of simple periodic orbits 
is performed and the largest Lyapunov exponent is numerically calculated. Section IIIII also contains the analysis of 
distributions of intersection points on the singular set. The existence of a weak noise observed in the Fourier analysis 
of intersections seems to be a complementary indicator of chaotic behavior. 



II. FRW MODELS WITH SCALAR FIELDS AS A SIMPLE INDEFINITE DYNAMICAL SYSTEMS 



The dynamical systems of cosmological origin have many special features which distinguish them from those met in 
classical mechanics. It is in fact the origin of some problems and controversy in understanding of chaos in cosmology 

mm. 

In this section we shall study the dynamical complexity of a simple inflationary models of the Universe, regarded 
as Hamiltonian dynamical systems. They appeared, for example, in Linde's chaotic inflation scenario, where inflation 
is driven by the vacuum energy of a single slowly rolling inflaton field |Tfij . 

The dynamics of cosmological models, allowing for an inflaton field, has been studied by many authors |!6lll7| . The 
idea of inflation which was introduced to solve some of the underlying problems in the standard big-bang cosmology 
becomes strictly connected with the existence of a scalar field which generates the period of an accelerated expansion 
of the Universe. In this case its energy density becomes dominated by the potential energy V((f>) of the scalar field 4> 
(the inflaton). Although the dynamics of inflation depends on the specifics of the models, the basic mechanism lies 
in the equation of motion which for a homogeneous scalar field (</> = 4>{t)) takes the form 

4> + 3Hj> + V'(cl>) + ^ct> = 0, (I) 
6 

where an over dot represents a derivative whit respect to the cosmological time t, and V = dV/dcj), R is the Ricci 
scalar here, and the last term vanishes for minimally coupled scalar fields (in general the last term is £,R<fi, where 
£ = 1/6 for the case of conformally coupled scalar fields). 

Here we deal with a single homogeneous and conformally coupled scalar field <j> with potential V{<f) on the FRW 
background. The same system was previously considered in terms of the original dynamical systems without using 
the conformal Jacobi metric by many authors in the context of inflation and chaos [18l ll9L |2fJ . Our goal in this 
paper is to point out different manifestations of the complex dynamics in terms of geodesies in the Jacobi metric. The 
motivation for such a study is to obtain an additional physical insight using the conformal metric and the behavior 
of periodic orbits. 

Our cosmological model assumes the FRW spatially geometry, that is, the line element is of the form 

d,s 2 = a 2 ( V ){-d V 2 + d X 2 + f 2 (x)(de 2 + sin 2 ddip 2 )}, (2) 
where < ip < 2tt, < 6 < it, < x < ^ an d V is the conformal time dt/a — dr/; a - the scale factor; 

{sinx, < x < f k = +1 
x, o< x <^fc = o 
sinhx, < x < 00 k = —1 

k is the curvature index. 

The gravitational dynamics is derived from the Einstein- Hilbcrt action 

S g = \Jd i x^{R-2K), (3) 

where g is the determinant of the metric; — g = a 4 / 2 (x) sin# and R is the curvature scalar for metric given by 

„ f d k 1 

where a dot denotes differentiation with respect to r/. We use the Misner-Thorne- Wheeler (MTW) convention [2l| as 
well as a natural system of units, such that h—c — 8irG — 1 and (27r) 2 = 1. In so far as Robertson- Walker symmetry 
holds, the scalar field should be homogeneous <f> = 4>(t). 
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W = t;{-{pI + ka 2 ) + (pi + k^ 2 ) + m 2 a 2 ^ 2 + A^ 4 + Aa 4 }, (6) 



The action for a conformally coupled (real) scalar field is given by 

S* = -\J d 4 xV=g{d^d^+±R0 2 - 2V{<f)} (5) 

where V{4>) — \m 2 cf) 2 + hcfft is the assumed form of potential for this scalar field. 

After integration over the spatial variables (J d 3 x — 2tt 2 is the conformal volume of an spatial hypersurface of 
constant curvature) and discarding total derivatives in the full action we obtain dynamical system with two degrees 
of freedom: a and a rescaled scalar field if) : <j> — > if) = yJTjQacj) with Hamiltonian 

2 1 - >/',.- ' r "" ~ '/'■• 

where m 2 is the mass of the scalar field, A is constant proportional by the factor of 1/3 to the cosmological constant. 

The evolution of the system should by considered on the Ti = energy level for vacuum cosmology or on the 
Ti = —p r fl energy surface if we add a radiation component to the energy-momentum tensor whose energy density 
scales like p r — p r ,oa~ 4 , where p r ,o = const p3 |. Let us note that system (j^J belongs to a larger class of dynamical 
systems which we call simple indefinite mechanical systems. 

By simple indefinite mechanical system we understand the triple (M, g,V), where M. is the configuration space 
carrying a metric g which defines the indefinite kinetic energy form K. = hg(u, u), u e T X A4, x e M.. V is the 
potential function V : M. — > TI. which is C°°, and g has the Lorentz signature (—,+,+,+). 

The "simple" in the above context means that dynamical system has the natural form of Lagrange function C = 
\gapq a q P — V(q), where a, (3 — 1,...,N and q and q are generalized coordinates and velocities respectively. The 
Hamilton function for such a system is of the form 

H(p, q) = \g al3 p aP p + V(q), Pa = g aS3 cf. (7) 

For the system under consideration (q 1 ,*? 2 ) = (a,ip). Because of our general relativity and cosmology application 
Ti = E = const g a f)q a ¥ 3 = 2(25 — V(q)). Therefore trajectories of the system in the tangent space of 1Z 2N with 

coordinates (q a ,q a ) are situated in the domain described by Q = {(q a ,q a ) G K 2N : g a m a q P = ||u|| 2 = 2(25 - V{q))}. 

In the tangent space T q {TZ N ) it is natural to distinguish three classes of vectors, namely, a vector u is time-like if 
||u|| 2 < 0, space-like if ||u|| 2 > and null if ||u|| 2 = 0. 

In the configuration space we distinguish three subsets 

V s = {q e K N : 25 - V(q) < 0}, 

V T = {q S K N : 25 - V(q) > 0}, (8) 
dV = {q e K N : E - V{q) = 0}. 

Note that set dT> is the boundary set because in its neighborhood we can always find points of T>g and T>t 

dV s = dV T = dV. (9) 

In the three distinguished domains the character of the vector tangent to a trajectory is determined by the Hamil- 
tonian constraint Ti = E. Therefore if a trajectory changes the domain, say T>s to T>t, it crosses dV and the tangent 
vector to that trajectory at the point q £ &D is situated on the cone determined by the kinetic form: g a pq a qP = E. 
In Ref. @ we can find the proof that in the case of B(IX) model trajectory crosses the boundary set V = 0. During 
each of oscillations that occur close to boundary set the solution is instantaneously Kasner and therefore V — an 
infinite number of times. 

The physical trajectories of the simple indefinite systems arc geodesies on pseudo-Riemannian manifold without 
the boundary (on which metric is degenerate) if we define the metric 

g a p = 2\E - V\g af3 

and reparameterize the time variable rj — > s [23| : 

^ = 2\E-V\. 
or) 

In our further considerations, we will demonstrate the advantages of analysis behavior of trajectories in a neighborhood 
of the boundary set dV of the region classically admissible for motion. It is also a surface of degeneration of the Jacobi 
metric. This surface and points of its intersections with the trajectories of the system contain interesting information 
about complex behavior. For simplicity of presentation of this idea we assume k — +1, i.e. vacuum and closed FRW 
model is considered. The domain (classically) admissible for motion, as well as the boundary set, are shown on Fig. 
miq 1 =a, q 2 = ^). 
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FIG. 1: The domain admissible for motion of the FRW closed system with a conformally coupled scalar field. In the shaded 
area trajectories behave locally unstable Kg(u, u) < where K is the Gauss curvature for the Jacobi metric. 

III. DIFFERENT EVIDENCES OF CHAOTIC BEHAVIOR 
A. Poincare sections 

It is clear that if a trajectory passes trough the boundary set dV 7 then tangent vector to the trajectory is null, i.e. 
it lies on the cone g a fiq a q^ — 0. The physical trajectories of the indefinite mechanical systems for given total energy 
E (zero for vacuum cosmology) are geodesies if we choose metric in the form g a p = 2(E ~ V)g a p in both domains T>s 
and T>t- On the boundary E — V the metric is degenerate, which is a source of obstacles if we define the property of 
sensitive dependence on initial conditions in terms of curvature invariants (the Gauss curvature in our case) . In both 
open regions T>s and T>t the Euler-Lagrange equation for the Lagrangian C = \g a pq a (f — V(q), ( ' = gj) assumes 
the form of geodesies equation after reparameterization of the time variable |2^| 

t->s:^ = 2\E-V\ (10) 

where s is the natural parameter defined along geodesies. 

The criterion of local instability of geodesies flow can be formulated as Kg(u, u) < 0, where u is the vector tangent 
to the trajectory and K is the Gauss curvature for the Jacobi metric. This domain is represented by shaded regions 
on Fig. ^ Due to the representation dynamics as a geodesies flow one can imagine fictitious free falling particle in 
both domains T>$ and T>t meeting the singularity (infinite curvature K) at the boundary set &D which play the role 
of a scattering surface. 

Fig. I2I3I4I shows three trajectories in the configuration space for three different initial conditions together with 
Lyapunov exponents in time, calculated in the standard way. Note that in all cases Lyapunov principal exponent 
goes to zero as r\ — > oo which corresponds to the singularity at the conformal time. This suggests that the system 
under consideration has no property of sensitive dependence on initial conditions which is characteristic rather for 
integrable systems. Another property which distinguishes the system from other systems of classical mechanics with 
conception of absolute time is that they do not posses property of topological transitivity which is crucial point for 
understanding of classical chaos conception pij . 

It is useful to choose degeneration line V(a, (j>) = as a Poincare surface. Let us concentrate on the simplest case 
of k = +1 and A = A = 0. Then <j) = ±a/vT + a? is an algebraic equation of the boundary set. On Fig. [5] one can 
see Poincare section on the surface V(a, 4>) — on the plane (a, <p) for different initial conditions. 
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FIG. 2: The trajectory in the configuration space for initial conditions ao = 0, <f)o = 4>o = 0.5 (do calculated from Hamiltonian 
constraint) and Lyapunov principal exponent which tends to zero which may suggest that the system is regular. 
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FIG. 3: The trajectory in the configuration space for initial conditions ao = 0, <f)o = <^o = 0.55 (a'o calculated from Hamiltonian 
constraint) and Lyapunov principal exponent. 



While in Fig. 0we find most trajectories as regular, there are, of course, chaotic trajectories. Some details of their 
concentrations in neighborhood of saddle points and forming stochastic layers is illustrated on Fig. 



B. Symbolic dynamics in detection of dynamical complexity 

Hadamard [2{| was the first to use methods of trajectories coding in investigations of geodesies on compact space 
with negative curvature, which belongs now to the field of symbolic dynamics. The significance of this method in the 
context of closed cosmology with a scalar field was pointed out by Kamenshchik et al. 26] . The crucial feature of this 
model is the existence of points of maximal expansion (d = 0, a < 0) and sometimes points of minimal contraction 
(d = 0, a > 0) or "bounces" . Then it is possible to classify all trajectories using localization of their points of maximal 
expansion and calculate the topological entropy which measures the growth of the number of periodical orbits as their 
period increases. Hence one can quantify the length of orbit by the number of symbols A - bounce of trajectory 
(d = 0, a > 0), B - crossing the line <j> = 0. Note that a > for physical reasons the extension of trajectories to a < 
domain means the extending the solutions beyond the big crunch which is only mathematically admissible |13|. 

For the model under consideration, two different coding procedures were used. In the first method we count all the 
intersections of trajectories with the axis {a = 0}. For cf) > we put symbol (1) while in opposite case if <j) < we put 



-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 5 10 15 20 25 30 35 40 

a t 

FIG. 4: The trajectory in the configuration space for initial conditions ao = 0, <f)o = 4>o = 0.65 (do calculated from Hamiltonian 
constraint) and Lyapunov principal exponent. 
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FIG. 5: Poincare section surface V(a, <f>) = and a < 0, (f> > 0. 



symbol (0). Another method of codding trajectories rests on the analysis of intersections points with the boundary 
set dT> defined as the surface V(a, <f) = in the configuration space ((/>, a). One can quantify the length of the orbit 
by two symbols : (1) - if <fi > at intersection point, and (0) in the opposite case if (f> < 0. In both approaches 
mentioned before the trajectories are represented by a sequence of zeros and ones. The next step in our analysis is 
the division of all coding trajectories into blocks in the simplest way. The blocks consist of to letters (00), (01), (10) 
and (11) and blocks (01) and (10) are treated as a identical. Then, after counting different blocks one can calculate 
the Shannon informative entropy following the rule 

ffs = y>lni (11) 



where Pi (i = 1, . . . , r) is the probability that a given block will appear in the trajectory coding. 
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FIG. 6: Details of Poincare section from Fig. 



The informative Shannon entropy characterizes the uncertainty degree of appearance of a given result. Following 
definition Ijllll H$ is a number which belongs to the interval [0,lnr]. If for all i pi assumes the same value equal 1/r 
then from definition we obtain Hs = hxr which is the maximal value of Hs- From Fig. [7|one can observe that 
Hs is a growing function of initial conditions <po to the limit which corresponds to purely chaotic behavior. 
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FIG. 7: Shannon informative entropy as a function of <f>o for two different codding methods. The horizontal line Hs = In 3 
denotes the limit for purely random process. 



C. Distribution of intersection points on the boundary set dV as a measure of complexity of dynamics 

It is interesting to check how the intersection points are distributed on the boundary line described by function 
4> = a/y/l + a 2 . Let L be the length along that line calculated form the origin to the intersection point and N denote 
number of such intersections in the interval L± AL. Of course ./V can be normalized to P and treated as a probability 
of finding a fictitious particle moving along a geodesic at a given point of dT>. The probability P as a function of L 
(normalized to unity) is shown in Fig. [S] For deeper analysis of the distribution of the points Fourier analysis was 
performed and the results of such analysis are shown in Fig. [5] The existence of weak noise in the power spectrum 
can indicate chaotic distribution of the intersection points. 

D. The existence of unstable periodic orbits as an indicator of complexity of dynamics 

It would by useful for classification of the periodic orbits to consider some reflection symmetries of the system. 
Of course, the system possesses reflection symmetry t — ► — t and q 1 — > — q 1 (q l — a, </>). Therefore it is sufficient to 
investigate motion in one quarter, say a > and <fi > 0, to reconstruct the motion in the remaining quarters. 

Due to the symmetries mentioned before the simplest periodic orbits in the configuration space can be grouped in 
the following five classes 

• trajectories of type 'I a' starting form the boundary line V(ao,4>o) = with initial conditions (drj,</>o) = (0,0) 
for which after a 1/4 of period the inflection point at ^ at o = is reached with (a, 4>) = (—0, 0) or (+4>, 0); see 

Fig. m 

• trajectories of type 'I b' also starting from the boundary set and reaching the singular point (a, <f)) = (0, 0) for 
which a = (f> or a = —(f) after a 1/4 of the full period; see Fig. ITT1 

• trajectories of type 'II a' starting from a point of configuration space (arj, 0) with (do, <f>o) — (0, ao) & n d arriving 
at the inflection point of cj>(a) diagram at a = (d, 0) = (—(f), 0) or (+4>, 0) after a 1/4 of full period; see Fig. ^| 
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FIG. 8: Intersection points of the boundary set 91?. Points of intersection of trajectory with set dT> as a function of time (left) 
and probability of finding a particle-universe on a given subset of the set dV during evolution (right) . 



Jk .kJL ,. .-J u 

50 100 150 200 250 ! 

frequency 



A 



50 400 450 500 50 



1 00 150 200 250 300 350 400 450 500 
frequency 



FIG. 9: Fourier analysis of distributions form Fig. [H| Existence of weak noise in the power spectrum can indicate chaotic 
distribution of intersection points. 
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Type 


go = ao 


Period 


I a 


1.90228463575087 
4.07891339562106 
6.14180160301955 


8.21179541505428 
9.68359976287486 
10.3389393707077 


I b 


2.63614545223571 
4.67384662462167 
6.68822113801299 


9.49181146895361 
10.4432886891063 
10.8883259280963 


II a 


3.08416398565973 
5.16022385714180 
7.19459890603674 


9.12958397138681 
10.06129848810491 
10.55370414489915 


II b 


1.64890985955010 
3.71764148794895 
5.72295574169836 


8.45628399932237 
10.07328719135456 
10.69967308122208 


III 


2.47587020635594 
2.75437747583008 
2.94236929747301 


17.99123630712047 
19.56534435565128 
20.15957644135879 



TABLE I: Unstable periodic trajectories. 
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FIG. 10: Periodic trajectories of type 'I a' in the configuration space. 



• trajectories of type 'II b' starting form the point (ao,0), (do,<Ao) = (0: a o) an d reaching the singular point 
(a, <f>) = (0, 0) and a — <fi or a = — (f> after a 1/4 of the full period; see Fig. ITS1 

• trajectories of type 'IIP (which are the union of all previously mentioned cases) starting from the boundary set 
V(ao, 4>o) — with (do, 4>o) = (0, 0), after a 1/4 of the full period they reach one inflection point at scale factor 
a, and ef> = with (d,0) = (0, ±a), after a 1/2 of the period they reach the symmetrical point to the initial 
condition with respect to 0-axis; see Fig. IbH 

Of course there are also non-symmetric periodic orbits present (see Fig. I15fl but they are not a subject of our 
consideration. In Table [I] we can find the periods of typical periodic orbits which we obtain by using modified 
multi shooting method 27]. To illustrate the property of sensitive dependence on initial conditions we investigate the 
evolution of the separation vector in the phase space in the term of "Lyapunov like" principal exponent. The principal 
Lyapunov exponent as well as the distance of the separation vector for initial conditions' separation Aa = 10 -12 , are 
illustrated in Fig. One can observe the existence of unstable periodic orbits in the model which should be treated 
as a strong evidence of chaos pgj . 
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FIG. 11: Periodic trajectories of type 'I b' in the configuration space. 
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FIG. 12: Periodic trajectories of type 'II a' in the configuration space. 



IV. CONCLUSIONS 

In the present paper the complex behavior of trajectories of FRW cosmological models with scalar field was inves- 
tigated by means of geodesies of the Jacobi metric. 

We pointed out the role of the singular set dT> of degeneration of the Jacobi metric in detection of complexity of 
dynamical behavior. We find that this domain can be useful as a Poincare surface. Moreover, the distribution of 
the intersection points as well as the existence of periodic orbits contain interesting information about the degree of 
complexity of its dynamics. Therefore all these investigations should be treated as a complementary description of 
chaotic behavior in a geometrical way. 

We have demonstrated the complexity of dynamics in the sense of (1) Poincare sections, (2) random distribution of 
intersection points, (3) the existence of unstable periodic orbits and (4) chaos in trajectories coding. All this evidence 
has rather mathematical sense because we prolong trajectories to the nonphysical domain a < 0. The true sense 
of complexity following Motter and Letelier 0] is in the nonintegrability of its dynamics. In Ref. [2^. I30L l3l| we 
have investigated the model under consideration in the framework of nonintegrability and we showed that they are 
non-integrable in the sense of nonexistence of meromorfic first integrals for the generic case of the model's parameters. 
To that end, the Yoshida, Ziglin, and strong methods of differential Galois group are used. 

The presented approach offers some new possibilities of coding trajectories for which a well defined conception of 
Kolmogorov complexity can be also applied |3^. In 1963-1965 A. N. Kolmogorov proposed to consider a measure of 
complexity in the framework of the general theory of algorithms. Let us consider model of dynamics in the form of 
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mathematical object that has binary string a as its complete description. Then we can use the length of the shortest 
program in bits as a measure of the complexity of the object |33j . 

Let us now consider a cosmological model as a mathematical object that has coding dynamics in the form of binary 
string a as its complete description. Then for the purpose of quantifying complexity, the notion of Kolmogorov idea 
can be useful. The question, whether the FRW cosmological model is complex in the Kolmogorov sense, is open. In 
our opinion the geometrical language introduced here offers a new, interesting possibility of trajectories coding in the 
form of binary strings which makes the answering this question easier. 
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FIG. 15: Non-symmetrical periodic trajectories in the configuration space. 
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FIG. 16: Stability analysis for the first trajectory from Fig. 1101 Principal Lyapunov exponent (left) and the distance of 
separation vector (right) for initial conditions' separation Aa = 10 -12 . 
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